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NON-UNIQUENESS OF ADMISSIBLE WEAK SOLUTIONS TO 
COMPRESSIBLE EULER SYSTEMS WITH SOURCE TERMS 

TIANWEN LUO, CHUNJING XIE, AND ZHOUPING XIN 


Abstract. We consider admissible weak solutions to the compressible Euler system 
with source terms, which include rotating shallow water system and the Euler system 
with damping as special examples. In the case of anti-symmetric sources such as ro¬ 
tations, for general piecewise Lipschitz initial densities and some suitably constructed 
initial momentum, we obtain infinitely many global admissible weak solutions. Fur¬ 
thermore, we construct a class of finite-states admissible weak solutions to the Euler 
system with anti-symmetric sources. Under the additional smallness assumption on 
the initial densities, we also obtain multiple global-in-time admissible weak solutions 
for more general sources including damping. The basic framework are based on the 
convex integration method developed by De Lellis and Szekelyhidi [Hill] for the Eu¬ 
ler system. One of the main ingredients of this paper is the construction of specified 
localized plane wave perturbations which are compatible with a given source term. 


1. Introduction 

The well-posedness of the compressible Euler system is one of the central issues in 
hyperbolic balance laws. The uniqueness of admissible solutions is of fundamental im¬ 
portance. Although there is a quite mature well-posedness theory for one dimensional 
hyperbolic conservation laws with small BV initial data the problem for multi¬ 
dimensional systems is very challenging. Recently, a major breakthrough for the unique¬ 
ness problem is made by De Lellis and Szekelyhidi in [T^fll] . Inspired by the surprising 
examples in [27] and [28], they developed a convex integration framework and obtained 
inhnitely many bounded weak solutions to the incompressible Euler system. The convex 
integration methods are later rehned to generate even Holder continuous solutions for 
the incompressible Euler system, see The ideas have also been applied to 

other systems of PDEs; see [9l|22l|29l[33] and the references therein. We refer to [15] for 
a general survey on the results in this direction. 

Multiple admissible solutions to the compressible Euler system are obtained in [U 
EHHilllilT], by an adaptation of the convex integration method. A major consequence 
is that the admissible weak solutions for the polytropic gases in multidimension are in 
general non-unique. These solutions are called “wild solutions” in the literature, which 
reflect the flexibility of the solution space with low regularity and are quite different in 

nature from those in one dimensional setting such as BH. It has been shown that many 
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of the available criteria Enn], with the exception of vanishing viscosity limit, are not 
able to single out a unique solution. 

It is interesting to investigate the stability mechanisms which may help to rule out 
the wild solutions. A natural candidate is the lower order dissipation and dispersion, 
such as damping and rotating forces. 

In this paper, we consider the compressible Euler system with source terms as follows 


dtp + V ■ (pu) = 0, 

dt{p\i) + V ■ (pu (g) u) + Vp(p) = B(pu), 


where {x,t) G x [0, oo) with or p, u, and p denote the density, velocity, 

and the pressure of the flows, respectively. Assume that the equation of states satishes 
p(0) = 0 and p'{p) > 0 for p > 0, and B is an it, x n constant matrix. In particular, 
the effects of damping and rotating forces are included in this model, where in the case 
of n = 2 and nondimensionalization, the matrix B has the form B = —I and B = J, 
respectively, with 


I 



and J 



( 2 ) 


Most of the previous investigations 0 El 0 [lain] focused on the isentropic Euler 
system corresponding to B = 0 in the system ([T]). The existence of inhnitely many 
bounded admissible solutions is hrst showed by De Lellis and Szekelyhidi in [Tl| for a 
special class of piecewise constant initial densities. The local-in-time existence of multiple 
admissible solutions for general smooth initial densities is proved in 0. When the initial 
density is close to a constant, the global existence of multiple admissible solutions is 
obtained in d- In these works 0[iaEi, the initial momentum are suitably constructed 
and not explicit. For some classical Riemann initial data connected by shocks, it is proved 
in 0171 show that the admissible weak solutions for polytropic gases in two-dimension 
are not unique. The non-uniqueness issue has also been studied in the class of dissipative 
weak solutions for the Euler-Fourier system [6] , the Euler-Korteweg-Poisson system [12] , 
and the Savage-Hutter model na, respectively, by adaptations of the convex integration 
method. 


The global well-posedness of classical solutions with damping and rotations has been 
well-known under certain smallness assumptions on the initial data; see for examples 
[3l[2ni[25l[35| . This is in sharp contrast with the Euler system whose classical solutions 
generally develop singularities no matter how smooth and small the initial data are; see 
for examples [26l[30|. The main aim of this paper is to investigate whether these lower 
order dissipations or dispersions can prevent the loss of uniqueness for admissible weak 
solutions. Our next aim is to construct wild solutions with special structures in order to 
gain a better understanding for the hne properties of weak solutions. 
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The weak solutions of the system ([T]) are the ones which satisfy ([T]) in the sense of 
distribution, i.e., 

Definition 1 (Weak solutions). A pair (p, m) G x [0, cxo); (0, cxo) x M”) is said 

to he a hounded weak solution of the system ([I]) with initial data 

p{x,0) = po and ni(x, 0) = ni'^(x), (3) 

if for any (<p,'0) ^ x [0, cxo);M x M"'), it holds that 



{pdtp + m ■ V(p)dxdt = — Po'p{x, 0)dx, 


rxi rxi 

m ■ dt'ip H-: + p(p) V ■ ■ Bm ) dxdt 

P 


m* • fj^x, 0)dx. 


It is well-known that weak solutions are in general not unique. Several admissibility 
criteria to exclude the non-physical solutions have been introduced. An important one 
is the so called entropy condition. Set 

= ( 4 ) 

Jo 

It is easy to see that (X(p) -f- (X(p) -|- +p(p))y) is an entropy-entropy flux pair 

for the system ([T]) [11]. Thus one can dehne the admissible weak solutions to the system 
([1]) as follows. 


Definition 2 (Admissible weak solutions). A hounded weak solution (p, n) is called an 
admissible weak solution to CO with initial data (13|) if for any non-negative test function 
(p G C^(R"' X [0, cxo)), it holds that 

[ [ (^ip) + + (lip) + + pip)] — • Vp + — ■ Bmcpdxdt 

Jo V V 2p y V 2p j p p 

+ [ (^ipoix)) + ifix, 0)dx > 0, 

Jm.” V ^Po / 

where X is defined in (j4)) . 


In fact, the inequality (]5|) is exactly the energy inequality for the system (JI]), and 
each strong solution to the Cauchy problem ([1]) and (jS]) satishes ([5|) with equality. Fur¬ 
thermore, if there is a strong solution for the problem ([1]) and ([3|), then any admissible 
weak solutions must coincide with it, see m- 

Now the main results in the paper can be states as follows. We start with the case 
in which B is an anti-symmetric matrix, such as in the rotating shallow water system. 
For piecewise constant initial densities, the following results hold. 
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Theorem 1. Suppose that B is an anti-symmetric constant matrix. Assume that = M” 
or T'^. Let po be a piecewise constant positive function in in the sense that there are a 
family of at most countably many mutually disjoint open sets Qi with the n-dimensional 
Hausdorff measure 'HA{VL \ (Ujfij)) = 0 and positive constants {pi] with 0 < infjpj < 
supj Pi < oo, such that 

Po{x)=pi forxeVti. 

Then there exists an m* G such that there are infinitely many global bounded 

admissible weak solutions (p, m) to the Cauchy problem ([T]) and ([3]). 

Furthermore, either of the following two cases holds: 

(1) (p, m)(x,t) = (pi,0) for a.e. {x,t) e Qi x [0, oo); 

(2) (p, m) has exactly Nf states in fli x [0, oo), where Nf = 

There exists at least one i such that the above second case holds in Qi x [0, oo). 

Several remarks are in order. 

Remark 1. As long as the piecewise constant functions satisfy the Rankine-Hugoniot 
conditions for the compressible Euler system, they are weak solutions for the system ([I]) 
associated with B = 0. However, if a piecewise constant function is a weak solution of 
the rotating shallow water system, it satisfies not only the Rankine-Hugoniot conditions, 
but also the algebraic system Bm = 0. It is guite surprising that TheoremU\ shows that 
we can still obtain families of non-trivial finite-states solutions in the class of bounded 
admissible solutions. This also implies that in the case B = J the finite-states weak 
solution (p, m) can not be continuous at any point in fli x [0, oo) where the second case 
in TheoremUl holds. 

Remark 2. Theorem{I\is also inspired by the study on the problem of finding deforma¬ 
tions with finitely many gradients in non-convex calculus of variations (see fM&f). 

Remark 3. For the class of Riemann initial data considered in 00'. we can also show 
the existence of infinitely many admissible finite-states solutions to ([I]) with B = 0. 

When the initial density is a general piecewise Lipschitz continuous function, the 
following results hold. 

Theorem 2. Suppose that B is an anti-symmetric constant matrix. Let po be a piecewise 
Lipschitz function in M”, then there exists an m^ G L°°(M"') such that there are infinitely 
many global bounded admissible weak solutions (p, m) to the Cauchy problem ([I]) and 
&■ 

Furthermore, there exists a T* > 0 such that each of these solutions (p, m) reduces 
to locally finite-states when t > T*, i.e. for a.e. {x,t) G M" x (T*,oo), there exists a 
neighborhood M of (x, t) such that (p, m) has exactly Nf states in J\f. 
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Remark 4. The densities of the solutions in Theorem\^in general develop discontinuities 
with complicated geometry even if they are smooth initially. Theorem\^ removes the key 
small-oscillation assumptions on the initial density in where the densities remain 
smooth for all time. 


When B is a general matrix, set 

/3 = max [0,sup{-e^Be : ^ e M", |e| = 1}] • (6) 

Our main results on the non-uniqueness of admissible weak solutions to the Euler system 
with general source term are as follows. 


Theorem 3. Suppose that 0 < p < p are two positive constants. There exists a positive 
constant e depending on p, p and (3. Given e G (0,£), if po satisfies 

0 < p< Pq< p 


and 


II (po - P**, Vpo)||l°°(T") < £ where p** = 


IT*^ 


Po{x)dx, 


( 7 ) 


then there exists an m'^ G such that the Cauchy problem ([I]) and Q has infinitely 

many global-in-time bounded admissible weak solutions which also satisfy 

\\{p- p\m)\\L^^jn) < (8) 

where k depends on e and tends to zero as e goes to zero. 


Remark 5. Theorem [3 is also true for the Cauchy problem in M"" provided that there 
exists a constant p such that po — p decays sufficiently fast at infinity. 

Remark 6. A weak solution to the Cauchy problem ([T]) and ([3]) is said to be dissipative 
if the total energy is non-increasing. Recall that admissible weak solutions satisfy the 
local energy inequality dS]), which implies the non-increasing of the total energy. The 
admissible criterion ([5]) is more restrictive than the dissipative criterion. In fact Theorem 
m holds without the smallness assumption o in the class of dissipative weak solutions. 
For the recent studies on uniqueness of dissipative weak solutions, see 


We now make some comments on the analysis in this paper. We adapt the convex 
integration framework in [HIH] and reformulate the system ([T]) as a linear system 
coupled with nonlinear pointwise constrains. Subsolutions are defined by relaxing the 
constrain sets. In order to obtain finite-states solutions we consider more general relaxed 
sets than those in [HIH]. The key for the convex integration scheme [131IM] is to analyze 
the wave cone Awhich generates localized plane waves of the linear system .if. Since 
the plane waves of are different when the source terms appear, the results for the 
Euler system with B = 0 in [T3],[T3] do not apply directly. Our key observation is that 
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in the high frequency regime, the plane wave solutions for the homogeneous system 
^0 are good approximations for . After correcting the errors by solving divergence 
equations [16], we find an elementary method to generate localized plane wave solutions 
to for any constant matrix B. In the case that B = 0, our method also gives an 
alternative construction for localized plane waves from the ones in [TSlIT^ . 

To obtain weak solutions, we use a convex integration scheme to iterate subsolutions. 
The scheme is inspired by [8l[3l| . One of the key ingredients of this paper is that we get 
some suitable relaxed sets with extreme points by a careful perturbation argument 

using the Caratheodory’s theorem for convex sets, which plays an important role in 
constructing finite-states weak solutions. When the initial densities have large variations, 
we construct non-smooth subsolution ansatz with complicated geometry building on the 
smooth ansatz in HT]. This allows us to remove the key small-oscillation assumption on 
initial densities in ini, and is the key to proving Theorem [2l The major idea to prove 
Theorem [3] is to construct a strict subsolution ansatz employing the acoustic potential 
T in [T7] . 

The rest of the paper is organized as follows. The reformulations of the problem 
and the concepts of subsolutions are introduced in Section O We analyze the localized 
plane waves in Section [3l In Section H] weak solutions are obtained by iterating the 
subsolutions using the localized plane waves. In Section [S] suitable subsolutions ansatz 
with various initial data are constructed and the main results are proved. 


Notations. The following notations are used in the rest of the paper. The bold letters 
in lower and upper cases are used to denote vectors and matrices, respectively. Set 
R+ = (0, cxo), R+ = [0, cxo), and S"'~^ to be the unit sphere of R". Denote by the 
vector space consisting of symmetric trace-free n x n matrices. Let 


Nn = dim(R” x 


n(n 3) ^ 

2 


and N* 


n{n + 3) 
2 


(9) 


Given a set /L C R” x denote by conv K the convex hull of K, and int conv K the 
interior of conv K. The space G^(S!);93) is the set of functions mapping from 2} to 93 
with continuous derivatives up to order k {k can be inhnity), and (7^(2); 92) is a subset 
of 92) for functions with compact support. Suppose / G G°(2); 92), let 


image(/) = {r G 92 : r = /(h) for some h G 29} (10) 

denote the image of /. A function / G L°°(R+ x R”) is said to be in G/oc(R-i-; -h))f*(R")), 
if for any test function 0 G L^(R") and any compact subset J C R+, 


fcjydx —)■ 


f(f>dx uniformly on J. 


( 11 ) 
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Denote by 

cm = {f\f e [ f = 0}. 

Jx) 

Let tr(A) denote the trace of the matrix A. Two matrices A and B are said to be 
A < B (or A < B) if for any ^ G S'" one has 

eA^<eB^ (ore^Ae<e^BO. 

Suppose that A, 5 G D which is a linear space, then [A, B] is called the line segment 
connecting A and B, i.e., 

[A,B] = {eA + {i-e)B\e e [o,i]}. 

Suppose M, N are two subsets in a vector space Y. Dehne 

M + N = {yeY:y = y' + y",y' G M,y" G N.} 

Given a space-time set D C then = {x\{x,t) G V} is the projection of V on 

R" and = V \ [Vt x {t}). Let ^ = (x, ^ be a point in R"+^ and denote by ^r{x) and 
Qr{z) = Qr{x,i) the open space cube and space-time cube, respectively, i.e. 

^r{x) = {x G R" : \Xi — Xi\ < -, i = 1, • • • ,77,} (12) 

and 

(5r(^) = { 2 : G R X R" : |t - t| < ^, |xj - Xi| < ^, i = 1, • • • , n}. (13) 

For simplicity, denote Qr = Qr(0). Denote by BA the fc-dimensional Hausdorff measure. 
Given a set O C R^, we also use \0\ to denote 1-A{0) when there is no ambiguity. 

2. Subsolutions 

In this section, the system ([T]) is reformulated as a differential inclusion in the spirit of 
De Lellis and Szekelyhidi [Hill] . Then subsolutions as relaxed inclusions are introduced. 
However, the formulation here allows for more general constrain sets compared with 
M, which will be employed to obtain hnite-states weak solutions. Our formulation is 
inspired by P IT01I33] . 

Given p > 0 and g > 0, set 

Kp^q = {(m, U) G R" X - U = ql}. 

By taking traces, it is easy to see that 

Kp,q = {(m, U) G R" X : |m|^ = npg, U = ™ ^ ™ — ql}. (14) 

Therefore, Kp^q is bounded in R" x We denote 

= sup |tc| < 00 . 

luE conv Kp^q 


( 15 ) 
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Observe that a bounded weak solution (p, m) to the system ([T]) is equivalent to a 
bounded quadruple (p, m, U, q) satisfying the differential equations 

dtp + V ■ m = 0, (16) 

dtm + V ■ U + V (p(p) + q) = Bm, (17) 

in the sense of distribution and the constrains 

(m,U)(a:,f) e 77p(a;,t),g(x,i) a.e. (18) 

The subsolutions are defined by relaxing the constrains as follows. 

Definition 3. The quadruple (p, m, U, q) G L°°(M"' x M+; x M” x x ]R_|_) is called 
a subsolution of the system m if it solves (uni) and 03 in the sense of distribution and 
satisfies the relaxed constrains 

{m,lJ){x,t) e conv a.e. m M” x M+. (19) 

Furthermore, given a family of compact sets K{x,t) C. Kp(x,t),q{x,t) such that the map 
{x, t) I—)■ K(^x,t) is continuous in a space-time open set V in the Hausdorff distance ( 
Definition 7.3.1.]), a subsolution (p, m, U, g) is said to be strict in D with constrain sets 
^(x,t) if 

(p,m,U,g) e CioM+-LZ{Rl). (p,m,U,g)l^ G C\V), 

and 

(m, U)(x, t) G int conv K(^x,t) inV, p > 0 and q>0inV. (20) 

Notice that the strict subsolutions defined in [11[T3] correspond to the cases that the 
constrain sets = Kp(x,t),q{x,t) in Definition [21 

The following lemma, appeared in [13] , plays an important role in constructing strict 
subsolutions, and implies in particular that (0, 0) G int conv for p, g > 0. 

Lemma 1. The convex hull of Kp^q is characterized as: 

conv Kp^q = {(n, V) G M" x :n®n-pV < pgl}. (21) 

Consequently, 

int conv Kpq = {(n, V) G M” x ■. n®n — pV < pgl}. (22) 

Proof. In the case of Jt'i i/„, the characterization fl2T]) is obtained in [T3|. The general 
case follows from the linear transform 

(m, U) ^ —U ) , for (m, U) G x (23) 

Vv/fW nq J 

which maps Kp g one-to-one to Adp/n, and preserves convex hulls. □ 
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In the rest of this section, we consider convex sets with hnite extreme points, which 
will be employed for constructing hnite-states weak solutions. First, the following ele¬ 
mentary result holds. 

Lemma 2. Suppose that K is a compact set in M®. Let C be a compact set such that 
C C int conv K. Then there exists a finite set K of K such that C C int conv K. 

Proof. For w & C, there exists a simplex S such that tc G int iS and S C int conv K. 
Let {iSj} be a hnite sequence of simplexes which cover the compact set C. It follows 
from Caratheodory’s theorem that each vertex Wj^- of Si can be represented as convex 
combinations of at most s-|-l extreme points {Vij^kYkLi of conv K. Since K is compact, 
the extreme points of conv K must belong to K. Let K be the collection of these extreme 
points Clearly that C C Ujint Si C int conv K. □ 

In the case that C consists of a single point, the following sharp result holds. 

Lemma 3. Let w G int conv Then there exists a finite set of Nf points fC C Kp^q 

such that w G int conv JC, where Nf = Nn -|- 1 = as defined in ([9]) . 

Proof. It follows from Lemma[2]that there exists a hnite set F = {wi}^^ C Kp^q such that 
w G int conv F. Using Caratheodory’s theorem yields that there is a subset 
with w G conv {wi^j^i- However, it may happen that w ^ int conv The 

strategy is to avoid the degenerate case by perturbing the vertices slightly. An Nn-set 
G = {vj}^Li C Kp^q is said to be non-degenerate if 

{vj — is a set of linearly independent vectors in x (24) 

The number of distinct A'n-subsets of F is J = Suppose that for some 

k with 0 < fc < J, the set C Kp^q has been obtained with the following 

two properties: 

(1) there are at least k non-degenerate A'n-subsets of F^; 

(2) it holds that w G int conv F^. 

For k = 0, the set Fq = F clearly satishes the above two properties. Suppose that we 
have obtained Fq, • • • , F^ which satisfy the above two properties. Let Gi, • • ■ , G; be all 
the non-degenerate A„-subsets of F^. Set Fi = F^ ii I > k. Otherwise the number of 
non-degenerate A^^-subsets of F^ is exactly k. Since k < (^), there is a degenerate Nn- 
subset G = of Ffc, i.e. the vectors {vi — linearly dependent. Consider 

a perturbation G = C Kp^q of G. It follows from a continuity argument that 

there exists a (5o > 0 such that w G int conv ((F^ \ G) U G) and that the A'n-subsets 
Gj = {Gj \ G) U G are non-degenerate for 1 < j < fc, whenever 

rnax \vi — Vi\ < Sq. (25) 
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In view of Lemma 0 ] below, it is not hard to see that there exists a set G which is non¬ 
degenerate and satishes fl2^ . Set ihc+i = \ G) UG. It is easy to see that F^+i C 

and satishes the two properties for k + 1. After at most (^) steps, we obtain F = Fj 
such that id G int conv F and all the A„-subsets of F are non-degenerate. 

It follows from Caratheodory’s theorem that there exists a A^*-subset )C of F such 
that w G conv /C. Furthermore, since any A„-subsets of F is non-degenerate, one has 
w G int conv /C. □ 

Lemma 4. Suppose p > 0 and q > 0. Given wq G Kp^q, for any given 5 > 0, the set 

Afs{wo) := {w G Kp^g : \w - WqI < 5} 

does not lie in any hyperplane V C M"" x 

Proof. First consider the case that p = l,g = 1/n, i.e. Kp^q = Ki^i/n- For a subset 
M cMP X §Q^"', let affdim(M) denote the dimension of its affine hull in R"' x i.e. 

the intersection of all affine sets containing M. For w G JG i/„, let 

^{w) = inf(affdim(A/ 5 (tc))). 

It is easy to see that Ki^i/n is diffeomorphic to via the map 

e ^ (G e e - -I) e for ^ g s ^-\ 

n 

Let Tu,(JG,i/n) denote the tangent space of iFi at w. It is easy to see that 

n-l = dim(T^(iFi_i/„)) < 2l(w;) < Nn, for w G 

We claim that equals to some positive constant Nq on Fi i/„, due to homogeneity. 
To see this, consider the action of SO{n) in R” x inspired by [8]. It is easy to see 
that for any R G SO{n), the linear transform on R"' x dehned by 

U) = (Rm, RUR') for (m, U) G R" x 

preserves affine hulls and transforms /G,i/n info itself. Then it holds that 2l(Rtc) = 2l(tc) 
for w G Ki^i/n- Suppose Wi and W 2 are any two points of Ki^i/n which are of the forms 

Wi = --I), G e = 1,2. 

n 

There exists a rotation R G SO{n) such that G = H follows that W 2 = PSwi and 
thus ^{ 102 ) = 2t(tci). 

We now show that Nq = Nn- Suppose Nq < Nn- Fix some wq G Ki^i/n- Then there 
exist a (5o > 0 and an affine subspace Lq of dimension Nq such that Afsoi^o) C Lq. Recall 
from Lemma [T] that (0,0) G int conv hence the affine dimension of /G,i/n is Nn, 

which is larger than dim(Lo) = Aq. It follows that the number 

6 = sup{(5 > 0 : J\fs{wo) C Lq} 
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is finite. From the definition of 6 and the compactness of 5A/j(iCo) C there exists 

a tci G dAfg{wo) with the property that J\fs{wi) is not contained in Lq for any <5 > 0. 
Since 2l(tci) = Nq, there exist a > 0 and an affine subspace Li of dimension Nq such 
that C Li. It is easy to see that Lq ^ Li and thus dim(Lo HLi) < Nq. Observe 

that for 5 < (5 which is sufficiently close to 5, the intersection Afs{wo) H Afsj^{wi) is a 
nonempty open set of Therefore, for w G J\fs{wo) fl A/ii(tCi) C Lq 0 Li, it holds 

that 2l(tc) < dim(Lo OLi) < Nq. This contradicts with the fact that 2l(-) is constant on 
It follows that Nq = Nn and the case for JO i/„ is proved. 

The general case for Kp g follows from the linear transform 0231) which maps Kp g one- 
to-one to and preserves hyperplanes. This finishes the proof of the lemma. □ 

3. Localized Plane Waves 

In this section, localized plane waves are constructed, which form the building blocks 
for the iteration scheme to obtain weak solutions. 

Dehne the linear operator 

^(n,V) = (V-n,atn +V-V-Bn)b (26) 

Observe that if (p, m, U, g) solves ffT6l) and ffT7|) and (n,V) satisfies 

if(n,V) = 0, (27) 

then (p, m -|- n, U -|- V, q) also solves ffT6D and (ITTD . The convex integration scheme will 
employ localized plane waves of . It needs to be shown that the wave cone Adefined 
below (according to [33]) is suitably large. 

Definition 4. The wave cone associated with L/f is a subset o/M”' x to which 

there exists a constant O > 0 such that for any (h, V) G there exists a sequence 
(nfc, Vfc) £ C'“((5i;M" x Sq^"') solving if) = 0 such that 

• dist{{nk,'Vk), [—(ft, V), (ft, V)]) —)■ 0 uniformly in Qi as k ^ oo, 

• (fik, Vfc) 0 m the sense of distribution as k ^ oo, 

• |(fffc, Vfc)p(a;, f)dxdt > OKft, V)p for all /c G N. 

The functions (ft^, W) are called localized plane waves of . 

Let =2fo denote the homogeneous operator corresponding to B = 0 and A^^ the 
associated wave cone. It is proved in [Tl| that if Wi ^ W 2 , for any Wi G Kp g with p > 0 
and g > 0, 

Wi-W2eA^g. (28) 

The proof of employs the key property that wi and W 2 is connected by a plane wave 
solution to A^^Q^n, V) = 0 as follows. 
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Lemma 5. Let A denote the subset of ML x such that for (n, V) G A there exists a 
{T,f) G M X such that (n, V)/i(rt + ^ • x) is a plane wave solution to =Sfo(n, V) = 0 
for any smooth function h. Then for w' and w" G Kp^g with w' ^ w" 

w — w G A. 


Proof. For (n',V'), (n",V") G Kp^g, set (n,V) = (n', V') - (n",V"). Let ^ G 
satisfy ^ ■ n = 0 so that f ■ fi' = f ■ fi". Denote r = —{l/p)f ■ n'. Then one has 


1 

rn + V • ^ = rn + -(n' ^ n' — n" (g) n") ■ f, = {t + -f ■ n'){n' — n") = 0. 
P P 


This hnishes the proof of the lemma. 


□ 


However, when the source B is non-zero, for wi,W 2 G Kp^g there may not be plane 
wave solutions to with profiles Wi — W 2 . It seems initially unclear whether the property 
fl28|) is still true for A_ 5 f in the presence of the source terms. Our key observation is that: 
in the high-frequency regime, localized plane waves with sources can be constructed as 
perturbations of the plane waves of ^o- This is stated as follows. 

Lemma 6 (Localized plane waves). Suppose that w,Wi,W 2 G M” x satisfy 

w = fiiWi + /X 2 IT 2 , Pi + P 2 = 1, Pi > 0, and w = W 2 — Wi E A. 

Given an open space-time set O and £ > 0, there exists a w E 0^(0; M” x So^”) 
satisfying the following properties 

(1) w{-,t) E £(M"') and w solves the eguation fl27D . i.e. 

= 0 (29) 

(2) dist{w -\- w{x, t), [iCi, W 2 \) < e for (x, f) G O; 

( 3 ) there exist two disjoint open sets Oi C O such that for i = 1,2, 

\w + w{x, t) — Wi\ < e for {x, t) E Op, < e. 

The following characterization of A_ 5 f is a consequence of Lemma El 

Corollary 1. A c A^. 

Proof of CorollaryUl Let w = (ti,V) G A. For each fc G N, applying Lemma El to 
ta = 0 = |(—ta) + with £ = ^ yields that there exists Vk E Cf°{Qi) satisfying fl2^ 
and 

f \vk\‘^dxdt -E- |tD 
Qi 


dist(Tfc, [——?■ 0, 


as fc —)■ 00 . 
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Let be the decomposition of Qi into mutually disjoint cubes of length 

2-K Let 

2 ^ 0 + 1 ) / \ 
ihfc = ^ af^ where af\z) = Vk e 

It is clear that ^Wk = 0 and dist(M)fc, [—'iv,w]) —)■ 0 as /c —)■ cxo. Since G £(R"') 

and diam(suppaj-^^) < \/{n + 1)2“^ —)■ 0, one has tCfc ^ 0 in distribution. Notice that 

2k{7i + l) 

[ iwfcp= Y1 f 

jQi J JQi 

Therefore, for sufficiently large k, one has 

j \wk\'^dxdt > -|tap. 

Jqi 2 

This shows that ra G and hnishes the proof of the corollary. □ 

The rest of the section is devoted to the proof of Lemma [6l 

To correct the errors from the source terms, we recall the linear operators which solve 
divergence equations in Sq^”) ( [16]) and adapt them to the whole spaces. 

Lemma 7. Proposition 5.1] There exists a linear operator TZ-^n from C°°(T";R"') 
to Sq^”) such that for any f G R”), it holds that 

V-77T.[f] = f-/ f, 

and 

ll’^T"[f]l|ci-“(T") < a)||f||c'“(T")- 

We also need the following lemma, whose proof is elementary and will be postponed 
to the end of the section. 

Lemma 8. For any f G C'))°(R"'; R”), there exists an 77[f] G C'°°(R”; Sq^”) satisfying 

V ■ 77[f] = f. 

Furthermore, TZ satisfies the following properties: 

(1) TZ is a linear operator from C)?°(R"';R”) to C°°(R"'; 

(2) 77[A^f] is a linear combination of third order derivatives o/f; 

(3) suppTZ[A‘^i] C suppi andTZ[A^f] G £(R"'), i.e., 
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(4) there exists a constant 0 < a < 1 such that 

||’^[f]||c“(R") < C'max(||f||il(Rn), ||f||£oo(Rn)), 

where the constants C and a depend only on the dimension n. 


and 


h{s) = 

|-h2, 




Proof of Lemma\^ Let 5 be a small positive constant to be chosen later. Fix a smooth 
cut-off function 0 G C^{0) such that 

(30) 

Set 

s e (0,/ii], 
s G (/xi, 1], 

and extend h as a periodic function with period 1. Notice that h{s)ds = 0. Let ho be 
a 1-periodic smooth approximation of h such that 

- P2 < ho < til, G [0,1] : h{s) ^ ho(s)}) < S, f ho{s)ds = 0. (31) 

Jo 

For k = 0,1, ■ ■ ■, let 

hk+i{s) = / hk{a)da, hk+i{s) = hk+i{s) - / hk+i{(T)da. 

Jo Jo 

Thus for any k > 0, one has 


hk\s) = hk-j{s) for 0 < j <k, I hk{s)ds = 0, 


and 


||hfc||L°° ^ ||hfc_i||Loo < • • • < llhollnoo. 

Denote w = (n, V) and let (r, .^) G M x be associated with tc G A. Set 
rt = X~^A^[nhe{XTt + Xf ■ x)(j)], V' = A“®A^[Vh6(Art -|- A^ • x)0], 


(32) 


and 


n" = -VA-^V ■ n', V" = 7^[B(n' + n") - dt{n' + n") - V ■ V'], 

where JZ and A“^u := AA * u are the operators dehned in the proof of Lemma [H and 
A G M is a large positive constant to be determined later. Dehne 


w = (n, V) = (n' + n", V' + V"). 

It will be verihed below that (12^ holds, suppta C supp0 and 

||tc(x, t) — who{XTt + Xf ■ x)(j){x, t)llioo(o) < C{w, (j), ho)X~^. (33) 

Indeed, it follows from Lemma [S] that 

V • n = V ■ (n' -I- n") = V • n' — V ■ VA“^V ■ n' = 0, 
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and 

V ■ V = V ■ V' + V ■ 7^[Br^ -dtU-V ■ V'] = Bn - dtu. 

Therefore, (n, V) solves the linear system =Sf(n, V) = 0. 

Direct computations show that 

n'= fihQ^Xrt + ■ x)(j) + \~^fi ^ Ca^pd'^ {hQ^Xrt + X^ ■ x)) 

l/3|>L|a+/5|=6 

= nho{XTt + • x)(f) + A“®n ^ (Art + A^ • x)d^4> 

l/5|>l,l«+/3|=6 

and 

V' = Vho(Art + Ae-x)0 + A-®V ^ A'^'Q^jC^e-iaKArt + A^ • x)af0. (35) 

l^l>l,l«+/3|=6 

It follows from fl32|) that supp(n', V') C supp0, and fj^„(n',V')dx = 0. The expressions 
(I3T|) and (l3^ imply that 

||(n', V') - (n, V)/io(Art + A^ •a;)(/)||Loo(o) < C'd^l, |V|,0,/io)A"d (36) 
It follows from the property n ■ ^ = 0 that 

n" = —A“®VA”^A^V ■ [nh6(Art + A^ • x)(j)] = —A“®VA^[(n ■ V(fyhQ^Xrt + A^ • x)]. 
Hence it holds that suppn" C supp0, n''dx = 0 and 

\\n''\\L--{o) < <^(1^1,0, ho)X~^. (37) 

Since B commutes with the Laplacian A, direct computations yield that 
V" = -n[dtn' + V ■ V' + dtn" - Bn] 

= —X~^TZ{A^[dt {fth^ (Art + A.^ • x) 0) + V ■ (Vhe (Art + A.^ • x) 0)] 

— cttVA^[(V0 • fi)h^ (Xrt + A^ • x)] 

— BA^[nh6(Art + A.^ • x)0] + BVA^[(V(/) • n)hQ{XTt + X^ ■ x)]} 

= -X-^n{A^[{ndt(j) + V • V0) he (Art + A^ • x)] 

— A'^V dt[{'V (j) ■ n)hQ (Xrt + A^ • x)] 

— A^B [A(nh6(Art + A.^ • x)0) — V((V0 • n)he{XTt + A.^ • x))]}, 
where rn + V • .^ = 0 has been used in the last equality. Hence V" can be written as 

V" = -x-^n[A^i], 

where 

f = A[(ncti0 + V ■ V0) he{XTt + X^ ■ x)] — Vctt[(V0 • n)hQ {Xrt + A^ • x)] 

— B{A[nh6(Art + A^ • x)0] — V[(V0 ■ fi)hQ{XTt + A^ • x)]}. 
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It follows from Lemma [S] that 7^[A^f] is a linear combination of third order derivatives 
of f. Therefore, suppV" C suppf C supp0, and \'''dx = 0. Combining the above 
computations yields fl2^ and that supptc C supp0. Furthermore, it holds that 

||V"|Uoc(o) = A-®||7^[A2f]|Uco(o) < C{\n\,\V\,<p,ho)X-\ 
which together with and fim) leads to ([23D- 

Since for {x,t) G 0,w{hQ(j)){x,t) G it follows from fl3^ that 

dist(t(; + w, [tci, wq\) < dist(t(; + whocj), [tci, W 2 ]) + |tc — who^Xrt + X^ ■ x)(f)\ 

< C{iu, 0, ho)A“^ 

Dehne the disjoint open sets Oi as 

= I (x, f) G : |tc + who^Xrt + A.^ • x)(j){x, t) — Wi\ < min 

It follows from (l30|) and (IM]) that can be arbitrary close to for 5 

sufficiently small and A sufficiently large. For {x,t) E Oi, it holds that 

\w + w{x, t) — Wi\ < ltd + who^Xrt + X^ ■ x)(j) — Wi \ + \w — who{XTt + X^ ■ x)(j)\ 

< I + C{w,(j), ho)A"h 

Therefore Properties (2) and (3) follow by hrst choosing 6 sufficiently small and then A 
sufficiently large. The proof of the proposition is hnished. □ 


/£ \W2 - Wil X 1 

V2’ 4 Jj- 


Proof of LemmalB Given f G dehne g = A where A = Af * u with 

the Newtonian potential Af dehned by 


Set 

^[f] 

where 


Af{x) 


1 1 


F(?7,) \ x \'^~‘^ 

— In Ixl 
27r ' ' 


for n > 3, 
for n = 2. 


^[vpg + (vpg)‘] + p (Vg)‘i - ^(V. g)i. 


Pg := g — VA ■ g. 

Then one can verify the properties of TZ by straightforward computations. First, one has 

Ti — 2 Ti 1 

V ■ 7^[f] = —--y[APg + VdivPg] + —-—[Ag + Vdivg]--(Vdivg) 

l{n — Ij l[n — Ij n — 1 

Ti — 2 Ti 1 

= 2 (TrT) Gs - Vdivg + 0] + + vdivg] - ^^^(Vdivg) 

= Ag = AA-^f = f. 
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It is clear that 7^[f] is a symmetric matrix and linear in f. Fnrthermore, 

= 0 + ^(V.g)-^(V.g) = 0. 

n — 1 n — 1 

Therefore, IZ maps to Sq^”). This proves the existence of IZ and 

the property (1). 

Since 

PAf = Af - Vdivf, 

one has 

71 — 2 

7^[A2f] = _ [V(Af - Vdivf) + (V(Af - Vdivf))*] 

^ (^33) 

This implies the second property of IZ. 

It follows from (j38|) that snpp7^[A^f] C snppf. If f e C“(M"';M"'), then taking 
integral for fl38l) on M"" yields 

[ n[A^i]dx = 0 . 

Hence one has the third property of TZ. 

Note that 

||f|Up(R'*) < max(||f||ii(]Rn), ||f||ioo(Rn)) for 1 < ^ < cx). 

Thns the fonrth property of 7Z is then an easy consequence of the boundedness of the 
Riesz operators and Sobolev embeddings [32] . The proof of the lemma is completed. □ 

4. Convex Integration 

In this section, weak solutions are obtained by iterating strict subsolutions, based on 
the method of convex integration. The main goal is to show the following proposition. 

Proposition 1. If (p, m, U, g) is a strict subsolution in T), then there exist infinitely 
many pairs (m, U) such that (p, m, U, q) are subsolutions and 

(m, U)(a:, t) G /or a.e. {x,t) inV, and suyvim — m, U — U) c P, (39) 

where K(^x,t) is the compact set associated with the strict subsolution (p, m, U, g) in Def¬ 
inition 0 
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To obtain weak solutions, we use a convex integration scheme to construct strongly 
convergent sequences of strict subsolutions. This gives a more constructive proof. The 
scheme is similar to the construction in [ 8 l| 23 ] , and the convergence relies on the method 
of controlled weak convergence in [3l] . An inductive argument using Lemma [6] yields 
the following lemma. 

Lemma 9. Suppose that K C Kp^q, and w G int conv K is a constant vector. For any 
£ > 0, there exists aw& x such that 

(1) w satisfies (I29j); 

(2) w + image{w) C int conv K, where image{w) = w{Qi) as defined in flTOll and 
image{w) denotes its closure in M"' x 

(3) the following integral estimate holds 

/ disfiw + w{z), K)dz < e. (40) 

jQi 

Proof. It follows from Lemma [2] that there exists a finite set C K such that 

w G int conv For fi G (0,1), dehne 

Lg = = w + {I - I3){wi - 

It is easy to see that for any 0 < /92 < /3i < 1, 

conv C int conv (Z int conv K. 

Observe that w G int conv Lg for any (3 G (0,1). Choose 5 such that 5 max* \wi — w\ < 
e/A. Set = Lq and 

^(i + l) ^ l^ii) y '■ w[ G W 2 — w[ E A, Ui E (0, 1), Z/l + 1^2 = 1}. 

Lemma [5] implies that 

= (1 — 5){wi — Wj) E A. 

It follows from Caratheodory’s theorem that 

Nr, 

= conv 
j=o 

where Nn = dim(M" x Sq^”). 

Let Tj = 2~^~^ for j = 0,1, 2, • • •. If tc G = Lq, set u) = 0. Clearly w E Lq C 
int conv Ls/ 2 - Observe that w = wf for some / and thus 

dist(tc, K) < max \wi — wf\ <6 max \wi — ti)| < e/A < (1 — To)b. 
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Assume inductively that for 0 < j < A;, if ta G and O C Qi is any open set, there 
exists awe C^{0) satisfying ([29]) and 


w + image(w) C int conv / dist(w + w{z), K)dz < (1 — Tj)e. (41) 

^ \^\ Jo 

If w e \ then there exist w[,W 2 e such that 

w = uiw'i + Z/ 2 W 2 , w '2 — a, z/j e (0,1), z/i + Z /2 = 1. 


L 


Given an open set O G Qi. Let eo he a small positive constant to be chosen later. In 


view of Lemma [H] there exists a. Wq E C^(0) satisfying fl29]) and two disjoint open sets 
Oi C O such that 


dist(w + Wo, [w(, W 2 ]) < Eo in O, 


and for i = 1, 2, 


w + Wo - w'l < £0 in Oi, < So- 


Since w' G C conv C int conv Lrf,s, it follows from the inductive assumption that 
there exists a Wj G C^iOi) satisfying (1291) and 



w[ + image (wj) C int conv 


(42) 


Let w = Wo + wi + W 2 . Clearly w G C^(0) and satishes ([29]). Furthermore, for 
ze0\(0iU02), 

w + w{z) = w + Wo('2) e i ?£0 + [w[, W 2 ] C + conv 
and for G Oi, i = 1, 2, 

w + w{z) = (w + Wo{z) — w[) + (w' + Wi{z)) G i?eo + conv 
Since conv C int conv Lr^.^^ 5 , for £0 sufficiently small it holds that 

Be^ + conv Lr^s C int conv Ln^^^s- 
Then w + image(w) C int conv 

Notice that ior z E Oi,w + w{z) = w + wo(^) + Wi{z). Thus for z E Oi, 


dist(w + w{z),K) < |w + wo(^) — + dist(w' + Wi{z), K) 

< eo + dist(w' + Wi{z), K). 


( 43 ) 
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Recalling the bounds flT^ and that Oi and O 2 are disjoint open subsets of O satisfying 

2 

i=l 

2 

i=l 

< 2eo, 


one can 


get that 

dist(t(; + w{z), K)dz < 


'o 


f dist{w + w{z), K)dz + 2C{p, q)\0 \ {UiOi)\ 

i=i 


< / 60 +dist{w[ +Wi{z),k)dz + dC{p,q)eQ 

i=l JOi 

— (I^ll + |C^2|)(1 ~ Tfc)e + CEq, 


where fj4^ is used in the last inequality. Since r^+i < r^, it holds that 
\ r _ (j£ 

W\ Jo ^ ^ < (1 - n+i)£, (44) 

for £0 sufficiently small. One can choose depending on {O), e, and r^+i so that 
it satishes fl43|) and (jHj). The induction step is completed. This hnishes the proof of the 
lemma. □ 


To deal with variable constrained sets K{^x,t)i one needs a stability result. 

Lemma 10. ( JT^ Lemma 1]) Let K be a compact set. For any compact set C C 
int conv K there exists a 5 > 0 depending on C and K such that for any compact set L 
with 

dH{K,L) < 6 (in Hausdorff distance) 

it holds that 

C C int conv L. 


In view of the above two lemmas, the following key perturbation property holds. 

Lemma 11. Suppose that {p, w,q) is a strict subsolution in a bounded open setF. Given 
e > d, there exists a compact set C G V and a sequence {w^} such that {p,Wk,q) are 
strict subsolutions in V and Wk — w G Cf°(V), suppiwk — w) G C,Wk ^ w in CL)f^. 
Furthermore, the following estimates hold 

/ dist{wk{x,t), K(x,t))dxdt < e. 

Jv 


( 46 ) 
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Proof. Set Co = 2 '?(^)) with €.{p,q) defined in (fT5D . Approximate V by 

a compact subset C' such that \ C) < For ( ^ V, applying Lemma |9] to 

ta(C) e int conv yields there exists v G Cf°{Qi) satisfying fl2^ and 

w{() + image(T) C int conv / dist(t(;(C) + v{z), K^)dz < —j^- (46) 

jQi 

It follows from Lemma fTOl and the continuity of w{z) and that there exists a r{C,) > 0 

such that 

w{z) + image(T) C int conv for z G Qr(c)(C)- (47) 

Furthermore, by shrinking r{C,) if necessary, the following estimates hold 

dH{K^,Ktf) < and |w(z) -'«;(C)| < ^ for z G Q^(C)(C)- (48) 

Let {C* = Qr(Ci)(0)}ili be a finite covering of the compact set C such that UC* C V, 
and let Vi G C“(Qi) be associated with Q. Set C := UCL Set tq = iminr((Cj). 

Apply the Whitney covering lemma |32] to the open sets Oi and {Oi \ (Uj=i ^j)}^ 2 ^ 
there exist a family of disjoint open cubes with 7f”+^(C \ (U^QO) = 0, each 

lying in some Por k = 1, 2, • • •, decomposing the cubes further if necessary, it 

follows that there exist finitely many disjointed open cubes {Q^ = Qrjizj)}jti, where 
Tj = Zj = depending on k, satisfying 

max T, < 2-*r„, «”+‘(C \ (Q Q’)) < 

j=i 

and for each j = 1, • • • , J*,, C for some i{j). 

Set aj{z) = G Cf°{Q^). Let Wk{z) = Y.jti «i(^) and set 


Jk 

Wk{z) = w{z) + Wk{z) = w{z) + ^ aj{z). 

i=i 

It is clear that Wk G Cf°{V) and supprcfc C UjQ^ C C. Since = 0, it holds that 
{p,Wk,q) solves flT^ and fITT)) . It follows from fl47)) and C that for z e , 

Wk{z) = w{z) + aj{z) G w(z) + image(Tj(j)) C int conv K^. 

For z E V \ (UjQ^), Wk(z) = w(z) G int conv K^. Therefore (p,Wk,q) is a strict subso¬ 
lution in P. 
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Substituting z = Zj + Vjz', it follows from and (HS]) that 


'Qi 


dist{wk{z), Kz)dz = / dist{w{z{z')) + aj{z{z')), K;^(^z'))dz' 


Qi 




dist(w(0(j)) + Vi(j){z'), + -^^dz' 


' Qi 
s IQ^' 


2\D\ 2\D 


■e. 


Thus it holds that 


'^k 


/ di'st{wk{z), Kz)dz = 2^ / di'si{wk{z),Kz)dz + / di'si{wk{z),Kz)dz 

I'D JVMvJiQi) 


D\{\JjQi) 


< + Col® \ (UjQ') I < I + Col® \ C'l 


< e. 


Let be the projection of on and denote by Zj = {xj,tj). Notice that 
ttj G C^iQ^), I aj{x,t)dx = 0, ||ail|L°° < ‘^Cq. 


(49) 


For 0 e it holds that 

Jk 

ijuk(pdx 


y / aj0dx 

,=i 


'^k 


/ aj{x,t){(f){x) — cl){xj))dx 


i=i 


< 2Co|supp0| max |0((r) — (l){xj)\. 


Since diam(,:S-^) < ^/n2 one has 


111 k(j) dx 


—)■ 0 uniformly in t as /c —)■ cx). 


It follows from density arguments that tCfc —?■ 0 in This completes the proof of 

the lemma. □ 


Proof of Proposition [3 Denote w = (m, U). Dehne 

= {w' G C/oc(lF+; : (p,w',q) are strict subsolutions in V and w' — w G C'“('D)}. 

Since w G X^, X^ is non-empty. It follows from flT5|) that X^ is a bounded set in 
L°°(EP X M+). Let X be the completion of X^ in Cioc{^+] L'2^,{R'^)) topology. Observe 
that X is metrizable due to the Banach-Alaoglu theorem and hence there is a metric dx 
such that (X, dx) forms a complete metric space (see [H]). It will be shown that the set 
of tc = (m, U) satisfying fl39|) is a dense set in X. 

Let be a family of bounded open subsets of V such that T>i <Z T >2 <Z ■ ■ ■ , and 

D = UDj. In the case that P is bounded, one can take Dj = "D for j = 1, 2, • • •. 
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Given w' G and Sq > 0. Set Wi = w' . Assume that for some k > Wk & 
has been obtained. Then {p,Wk,q) is a strict subsolution in Vk+i and thus Wk{z) G 
int conv Kzioi z E D. It follows that dist(tafc(z), Kz)dz > 0. Applying Lemma ITT] to 
Wk with Pfc+i, there exists a sequence {uj} satisfying that (p, Uj, q) are strict subsolutions 
in 'Dk+ii with 

Vi-WkE C^(Vk+i), Vi Wk in as i oo, 


and obey the estimates 


''V>k+i 


dist(Tj( 2 ;), Kz)dz < min <2 , - / dist(tCfe(; 2 ), Kz)dz 


IVi 


(50) 


Since Vi{z) G int conv for 2; G Vk+i and supp(Ti — Wk) C Vk+i, it follows that 

Vi{z) E int conv for z eV. It is easy to see that {p,Vi,q) are strict subsolutions in 

V. Since Vi — w = {vi — Wk) + {wk — w) E C^{V), it holds that Vi E X^. Since Vi Wk 

in as i —)■ cx), it is easy to see that for i sufficiently large, 

2 ' 


max 

l<j<k 


{vi - Wk)wkdz 


'Di 


< min < 2 


-k 


'Vi 


dist( 10 ^( 2 ;), Kz)dz 


(51) 


’ IOOIPjI 

Recalling that dx is induced by the topology of for i sufficiently large it holds 

that 

dx{vi, Wk) <2-’^eo. (52) 

Set Wk+i = Vi for some large i such that Wk+i satisfies flSTl) and fl52|) . Thus we obtain a 
sequence {wk} C X^ satisfying fl50|) . (I5TD . and fl52l) with n, replaced by Wk+i for k > 1. 

It follows from (l5^ that {wk} is a Cauchy sequence in X and thus Wk converges to 
some w E X. We now show that {wk} converges strongly in L‘^{Dj), for any hxed j. 

Fix any j > 1. For k > j, since Vi C Vj C Vk+i, it follows from (15(11) that 


\wk+i-Wk\dz> / dist(w;fc, - / dist(w;fc-M, 


'D, 


'D, 


'Di 


> 


dist(tCfc, Kz)dz — 


dist{wk+iiK^)dz 


'D, 


' D 


fc+1 




dist(tCfc, Kz)dz. 


Di 


Thus applying Holder’s inequality yields 

1 




\n 


\wk+i-Wk\dz\ > 


'Di 




\wk+i - Wk\dz 


'Di 


This together with (IHT]) . yields that 

WWk+lWh^Dp - W^kWh^Dp = \\Wk+l - Wk\\l 2 ^D, 


{wk - Wk+i)wkdz > 0. 


'Di 







24 


NON-UNIQUENESS OF EULER SYSTEM WITH SOURCE TERM 


Thus is a non-decreasing sequence. Since ||tCfc||L 2 (Dj) < C{p,q)\Vj\ is 

uniformly bounded, is a Cauchy sequence. It follows from fIFID that for 

k > m > j, the following estimates hold 


\\Wk — Wm\\‘i2(^D.) < 2{\\Wk — ITfc_i||^2(£,^.) + • • • + llWm-hl “ 1 1l2 )) 


k-1 


= 2 ^ \\wi+i\\l 2 ^jj.) - ||wz||i 2 (^.) - 2 / {wi+i - wi)widz 

l=m V J 

< ‘^{\\Wk\\L2{Dj) ~ \\'^m\\L2(Dj)) + ^ ■ 

Hence Wk converges strongly in L‘^{T>j). Then one can conclude from (j50l) and (j52i) that 
w(z) G Kz a.e. in V and dxiw^w') < Sq. 


Since the scheme holds for any w' G and £ > 0, the set of tc = (m, U) satisfying (1391) 
is a dense set in X. Since w ^ X^, {wk} is an infinite sequence in and thus a dense 
set in X is also inhnite. The proof of the proposition is hnished. □ 


As in [UlIlllIT], it is possible to construct weak solutions with initial data satisfying 
(l56|) . which are useful for constructing admissible solutions. It is straightforward to 
extend Definition [3] to the cases where the time domain M_|_ is replaced by [—r, cxo) for 
some r > 0. 

Corollary 2. Suppose that (p, m, U, g) is a strict subsolution m D C M” x [—r, cxo) for 


some T > 0, such that Dt=o C M” zs a non-empty open set and p{-, 0) = po with 


m ^(x, 0) = npo{x)q{x, 0), 

for a.e. x G (M"" \ T>t=o), 

(63) 

(rn, U)(x, f) G K p(x,t),q{x,t)i 

for a.e. {x,t) G (M” x M+ \T>). 

(54) 

Then there exists an m* G L°°{W^) with 

m'^(x) G (A'(a;,o)|R"), 

for a.e. x e A=o, 

(55) 

and m*p(a;) = npo(a;)g(a;, 0), 

for a.e. x E M", 

(66) 

such that the Cauchy problem ([T]) and ([3]) 

has infinitely many weak solutions (p, 

m^) 

satisfying 

m\x,t) G (A'(^,t)|R"), 

for a.e. {x, t) G V, 

(67) 


and \m!'\‘^{x,t) = np{x,t)q{x.t), for a.e. (x, t) G M"'x M_,_, (58) 

where iC(x,q|R'* = {n G M"' : {'n,\) G K(^x,t) for some V G 

Furthermore, there exists a sequence of divergence-free vector fields {liij} such that 

ihj G Cf°{V) and m-|- irij —)■ m*' in Czoc(K+; (59) 
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Proof. Denote w = (m, U). Define 

= {w G Cioc{[—T, oo); : (p, w, q) are strict subsolutions in V and w—ui G 

Given w G > 0, and a bounded open subset Dq C 'Dt=o- Modifying slightly the 

proof of Lemma [TTl one can obtain a sequence of functions {tCfc} C X^, which satishes 
dx{wk, w) ^ 0 and 

supp(tCfc — tc) C Do X / dist{wk{x,0), K(^x,o))dxdt < e. (60) 

J Qq 


For G Do, it follows from Lemma 0 Lemma [TUI and the continuity of w{z) and that 
there exist v G C^{Qi) and r(^) > 0 such that v satishes (12^ and 


/ dist{w{f,0) + v{x,t), K^^o) dxdt < 

'qi 4|iZo| 


Furthermore, for z G Qr(c)(^,0), we have 


dniKz-, K{i,o)) < 


w{z) + image(T) C int conv 
w{z) - w;(^,0)| < 


Since Qi = 


i Pn 
2’ 2 J 


X 


8|Do|’ 

, |], it follows from (iQTjl that 


8|Dc 


1 1 
' 2 ’ 2 


(61) 


-- -1 
2 ’ 2 J 


dist(w;(^, 0) + v{x, s), iF(^,,o))da; < 


4|Do|’ 


for some s G Let C be a compact subset of Do with 'H'^{Vt=o \ C) < where 

Go = 2suppG(p, g) with C{p,q) dehned in flTHl) . Denote by the space cube 

of length r(^) centered at Let {O* = be a hnite covering of C with 

Vi G C^{Qi) and Si associated with Set vq = minjr(^j). For k E N, there exists 
= .^pk){xf^)}jti such that \ (lj/=i ■^^)) < ^ 


< min(2 2 ^e), C for some i(j), for j = 1, • • • , 


Set 


i(x,t) = Vi(j) 


X - x^^^ t 

^ ■’ XP) + I ^ C!^i'^rj{Xj) X (-£,£)) 


r 


(k) 


j 


and let Wk = w + is easy to verify that Wk E and that (|60l) holds as in 

the proof of Lemma [TT] 

Using (1601) instead of Lemma ITTl the proof of Proposition [T] can be adapted to obtain 
{y;(^)} (y with vjd^ = (m, U), satisfying the estimates 


^[k+i) ^ ^(k) ^ ^ ^_2-k^ 2 -^), < 2 -^ 
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and that for any compact subset C C Vt=o^ 0)} is a strongly convergent sequence 

in L‘^{C) and 

dist(tc^*^^(x, 0), K(^x,o))dx ^ 0 as fc ^ oo. 

Let w' = (m', U') ;= linifctt;’^^) G X. Set m*(a;) = m'(a;, 0). It is easy to see that 

(m', U')(a:, 0) G -ft"(x,o) T^t=o and supp(m' — m, U' — U) C V. 

Since w' = for t ^ (—2“^, 2“^) for any A; > 1, {p,w',q) is a strict subsolution in 
Applying Proposition [T] to {p,w',q) in yields infinitely many pairs 

(m^U') G X, (m^U^)(•,0) = (m', U')(-, 0), (m', U')(a:, f) G a.e. in V. (62) 

It is clear that (l55l) and (l57)) are satished. Since K{x,t) C Kp{x,t),q(x,t), it is easy to see 
that (1561) and (1581) follow from (l53l) and (l54|) . 

Since (m'’,U'’) G X, there exist (mj,Uj) G X^ converging to (m^U^) in X. Let 
liij = nij — m. It follows from the dehnition of X^ that one has rhj G C^{V). Since 
rrij and m both solve (IT^ . irij are divergence-free vector fields. This finishes the proof 
of the corollary. □ 



5. Proof of the Main Results 

5.1. Constructions of finite-states admissible solutions. The starting points of the 
constructions are the piecewise constant stationary strict subsolutions inspired by ng. 

Proof of Theorem [3 Let y be a constant satisfying that y > maxjp(/9j). Let 

Qi = -P{pi) + X and I = {i : ft > 0}. 

For each i ^ I, Lemma [3] yields a set of Nf states /C* C Kp.^g. such that (0,0) G 

int conv /Cj. For i ^ I, let /Cj = {(0, 0)}. Set 

q{x,t) = fi, and K{x,t) = Xi, for x G hlj. (63) 

Let V = (U-^^ Rj) X (—l,cx)). It is easy to see that (p, 0,0,g) is a strict subsolution in 

T> with the piecewise constant constrained set K(^x,t)- Corollary [2] yields infinitely many 
weak solutions (p, m^) to the system ([T]) for some initial data (po, m^). 

For i ^ J, g = 0 in Rj. Hence it follows from (156|) and (|58|) that 
(p, m'')(a;, f) = (pi,0) for a.e. (x, t) G Rj x ]R_|_. 

For i G /, denote /C* = and w = (m'’,-^-gl). It follows from (162]) that 

TV* — ^ 

w G X and w{x, t) G {ih a.e. in Rj x M+. Hence w has at most Nf states in Hj. Since 
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/ N* 1 N* 

0 G int conv one can uniquely represent 0 as a convex combination of 


N* N* 

-''n ■‘'n 

0 = '^fiiw\ 0 </xz < 1, y^/iz = 1. (64) 

1=1 1=1 

It follows from the proof of Corollary |2] that there exist {wa:} C with 

'^k 

Wk^w in CL^^, wi = 0, Wk+i -Wk = '^ for k>l, 

i=i 

where satishes suppa^ C C x M+ for some l{j) and J aj{x,t)dx = 0 due 

to fl49|) . Hence Wk(x,t)dx = 0. For any compact time interval J C M+, it holds that 


w{x, t)dxdt = \im. / Wk{x,t)dxdt = 0. 


' QiXj 


' QiXj 


1 N'" 

Since w{x,t) G {tD'};d\a.e. in x M+, it holds that 


0 = / w{x,t)dxdt = where ui = ({(x, f) G x J : w{x,t) = w''}) . 

JUiXj 

It follows from fl64|) that z/; = x J| > 0. Hence (p, m^) has exactly N* states in Hj. 
It remains to consider the energy inequality ([5]). It follows from (]56|) and (15^ that 

lo Iq ^ ^ /j ^ 

Furthermore denote P{p,q) = p~^(X{p) + ^q + p{p))- It holds that 

r / Im^P \ /‘°° f 

/ ( X(p) + -h p(p) ) — • Vifdxdt = / P(p, g)m'’ • Vpdxdt. 

Jn\ 2p J p Jo 

It follows from (l5^ that there exist divergence-free vector helds liij G CJ°{TP) with 
liij —?■ in From the dehnition of P, it holds that rrij = 0 on dVLi for i G / and 

iiij = 0 in Hz for I ^ I. Hence one has 

rOQ n _ nOO n 

/ / F(p, q)mj ■ Vpdxdt = V' P(Pi, qi) / / • Vpdxdt = 0. 

Jo Ju jgj Jo J Qi 



Therefore, 




, g)m*' • Vpdxdt = 


lim 

3 




, q)va.j ■ Vpdxdt = 0. 


Since B is antisymmetric, it follows that (p, m'’) satishes ([5]) with equality. 


□ 
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5.2. Constructions of admissible solutions with a general source. Given po sat¬ 
isfying ([7]), the acoustic potential in m is employed to construct a strict subsolution 
such that p becomes constant in hnite time. 


Proof of Theorem\^ Consider hrst the case: po = P^- Let q = e for some ho > 0, 
where {3 is dehned in (jb]). It follows from (|22|) that (p**,0,0,g) is a strict subsolution 


Since V ■ = 


in P = X (— l,cxo) with constrain sets Kpt, q(^^ ^y Then Corollary [2] yields inhnitely 
many weak solutions (p**, m*’) to the system ([1]) with initial data (p**, 

—dtp'^ = 0, it follows from (El), fl56l) . and fl58l) that 


d, {!(/) + 


m 


b|2 


pB 


+ V- 


I(p’) + 


1 |m 


b|2 


pB 


m 


pB 


-B 


m 


• m 


n , 


pB 


<-{dtq + 2^q) = 0. 


Therefore, the energy inequality (ED holds for the weak solution (p**,m'’). 

Consider now the case that po is a small perturbation of p^ under the assumption ([7D. 
Let h{t) be a smooth cut-off function satisfying that 


h e 1)), h(0) = 1, 0 < h{t) < 1, and \h'{t)\ < 2. (65) 

Dehne 

p{x,t) = [l-h{t)]p^+ h{t)po{x), q{x,t)=p(p^)-p{p{x,t)) + x{t), (66) 

and 

m(a:,f) = h'{t)V^{x), U = 7^T-[-h"(t)VT(x) + h'(t)BVT(x)], (67) 

where x(t) is a function to be determined and T is the unique solution in /^(T'^) of the 
problem 

AT = p** — Po in T”. 


It follows standard estimates for the Poisson equation [12] and Lemma [7] that 


||d'||c 2 .“(T") A Ce and ||U||c'i,a('f>^) < C{\h'{t)\ + \h"{t)\)e. 


( 68 ) 


Furthermore, the density can be estimated as 

0 < p < p(a:,t) < p. Up** - p||l°°(T") < h(t) sup |Vpo| < h(t)£:. (69) 

The direct computation shows that (p, m, U, g) solves ffTBl) and flTTD . One has 

- U - gl = -^VT 0 VT - U - (p(p“) - p(p) + x) L 

It follows from fEHP and fEop that the inequality fl22l) holds, provided one can show 

X{t) > C{p,p)[hftfe+\hft)\ + \h%t)\ + h{t)]e (70) 

for t G (— ro,cxo) with some small tq > 0. Then (p, m, U, g) is a strict subsolution in 
P = T” X (—To, oo) with constrain sets Kp(^x,t),q(x,t)- Therefore Corollary [2] yields inhnitely 
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many weak solutions (p, m'’) to the system ([T]) with initial data It remains to 

check that these weak solutions are admissible. 

Recalling that {3 is defined in ([B]), it follows from flSB]) and fISS]) that 


+ 


m 


b|2 


2p 


+ V- 


X(p) + 


m 


b|2 


2p 


+ p(p) 


m 


P 


b" 


m" ■ fBm^ 


P 


- dtp{p)] + dtlip) +n(3[p{p^) - p{p{x,t)) + x{t)] 


V- 


'Ap)+P{p) + {n/2){p{p^) -p(p)). 


-m 


P 


n 


f 


+ ^xWV.^- 


(71) 


For t G [0,cxo), since supph C (—1,1), it follows that (p, m, U, g) = (pl*,0,0,x) for 
t > 1, thus reducing to the case of constant density. Hence if x(l) > 0, then x(t) can be 
defined on [1, cx)) satisfying flTOj) and the energy inequality (jS]). 

For t G [0,1), in view of (1^ . (1^ and the estimates 

|V-m‘'| = \dtp\ < \h\t)\e, |p(p“)-p(p)| <C'(p)(p“-p) < ^(p)^, 

|m^| = ^/np[p{pi) -p(p) + x] < C'(p)(v^ + ^/x), 


it is not hard to see that the weak solution (p, m'’) satisfies 


(X(p) + 


m 


b|2 


2p 


V- 


Ap) 


m 


b|2 


+ p(p) 


m 


m*' ■ (Bm^) 
P 


2 p ■ ' y p \ 

as long as y satisfies 

dtxit) < + x + x'''^+ '^)c forie [0,1). 


< 0 , 


(72) 

(73) 


Therefore there exists a e depending on p, p, and (3 such that if e < e, it is possible to 
find a smooth function y(t) satisfying (1701) and fl73l) on [0,1]. Hence the solution (p, m'’) 
satisfies the energy inequality (|5]). Observe that for sufficiently small e, one may take 
6 = y(0) to be small as well. Finally, the estimate (IH]) follows from the constructions. □ 


5.3. Constructions of admissible solutions with general density. The aim is to 

remove the small-oscillation assumption ([7]) on po in the case that B is anti-symmetric. 
However, since (1731) is a Riccati-type ordinary differential inequality, the assumption 
seems essential for the smooth subsolution ansatz in im. The key idea is the construc¬ 
tion of a non-smooth subsolution ansatz based on the ansatz in nzi, which transits to 
piecewise constant states before the blow-up time. 

Proof of Theorem [B It is assumed that po is piecewise Lipschitz. Therefore there exists 
a family of mutually disjoint open sets C M"' with po G 

oo / oo 

M^ = ljn; and -H" IJda 

i=l Vi=l 



0. 


(74) 
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and there are positive constants p, p, g snch that 

0 < p < inf po{x) < sup po{x) < p < oo and sup || Vpo||L°o(op <'g < oo. (75) 

i 

Let T G (0,1) and i? G (0,T) be two positive constants to be determined. It follows 
from the Whitney covering lemma [32l Theorem 3, p.l6] that there exits a countable 
family of mutually disjoint open cubes ( 2 ;(i))}^i g^ch that 

OO OO 

c ^i(j) for some i(j), 
i=i 


2=1 


and 


Tj < min 


dist(^^'’\ (1 + "2 (1 + ^?) 

Since every lies in some fli(j), so po G 

Let Tj be the unique solution in to the problem 

ATi = pf-po in 


for j G Id. 


du 


= 0 


on 




where 


Pi = 






Po{x)dx. 


Let h be a smooth function satisfying fl65|) and dehne 

hT{t) = h{t/T)eCr{{-T,T)). 

For {x,t) G X (—l,cxo), set 

Pi{x,t) = [1 - hT{t)]pl + hT(t)po{x), ini{x,t) = h'rp(t)V^i{x), 

and 


( 76 ) 


qi{x, t) = -p(pi(x, ()) + xit), 

where x(t) is a positive function to be determined satisfying that 

X{t) = X{T) > p(p) for t G [T, cx)). (77) 

Recall that 'H"(M"' \ (Uj^^*^)) = 0, it suffices to dehne p, m, q on each cube as 

p = Pi, m = nii, q = qi, in x (-1, cx). 

Let $(p) = i{x'^'‘"> + Tiy) and f{y) = p\ — Po(3^*'*^ + ^iV)- Then d) solves 

A<F = / in 


( 78 ) 
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Since 

||/llci(^i) < ^^illVpollL-(^W) < (1 + 

it follows from the standard estimates for the Poisson equation ng that 

Hence it holds that 

llV^.IUoo(^W) < r,|lV<h|Uoc(^,) < ^7(1 + 

Thus 

||ni||ioc(Kn) = |hy(f)| sup ||VTi||^oc(^(i)) < ||94m||ioc(]R™) < CT~^'d^. 

i 

Since the cubes {■^^} have disjoint interiors, so 


(80) 


j 

oo 


<J2 E iivuiU“i^“i+ E iivLiU“i^“i 

1=0 2*<|xU)|<2*+i |xU)|<l 

oo /oo 

<^(1 + | 2 f )-(”+^)i? 2 |H 2 *+i+i| +i 9 ^|H 2 | < (Pi^M ^ 2 “' + 1 I < C'-dl 


1=0 


1=0 


Thus 


m 


\h'T.{t)\ E II V'l'illi.ij.M) = CT-'tfl ||am|U.|,„, < CT-V 


Hence we can define 

U = 7?.[—^iin + Bm], 

It follows from Property (4) of Lemma [H] that 

||U||c7.(K.) < CT-^§\ (81) 

where the constant C depends only on a,n and the constant matrix B. Since 

z/■ mi = h' 2 .(f)z/• VTi = 0 ona^®, (82) 

and that dtpi + V ■ m* = 0 inside it holds that for p e C'^(M"' x (—1, cxo)) 


/ OO n _ POO p 

/ pdtp + m ■ V(p dxdt = / / 

-1 Jr” ■ J-1 J MO 


Pidt^p + mj ■ Vpdxdt = 0. 


Furthermore, direct computations yield 

dtm + V- U + V(p + g) — Bm = dtm + V ■ TZ[—dtm. + Bm] + Vx(f) — Bm = 0, 


where p + g = x is used. It follows that (p, m, U, q) solves flT^ and f[T7|) in the sense of 
distribution. 
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Since supp/i^ C (—T, T), denoting = x{T) — p(Pi), it holds that 

(m, U) = (0, 0) for t G [T, cx)), (p, q) = (pf, gj) for (x, t) G x [T, ex). (83) 

Since ™ ^ — U — gl < + |XJ| +p(p) — x(^)^ h follows from flHOj) and flHTj) 

that the inequality fl22]) holds provided that 

x{t) > p(p) + C{p, p)T-^i^^l(^_T,T){t), (84) 


where l(_T,r)(^) is the characteristic function of (—T,T). Suppose that fIMD holds in 
(—ro,T) for some small Tq > 0. Let 

P = X (-ro,T), X (T,x), V = Vul\JV^^'> 

j \ i 



For each i, Lemma [3] yields a set of N* states /C, C Kj j such that (0, 0) G int conv fCi. 

’’"i 

Set 

[iCi in F>W. 

It follows from flH^ and fl841) that (p, m, U, g) is a strict subsolution in "D with constrain 
sets K{x,t). Therefore Corollary [2] yields inhnitely many weak solutions (p, m*’) to the 
system ([T]) with initial data (po,m*). 

For f > T, in view of (l83D . as in the proof of Theorem [H we can show that (p, m) has 
exactly N* states in x [T, x), and that the energy inequality ([5]) holds with equality 
in [T, x). 

For t G [0,T), observe that in view of fl5^ and ([82]), it suffices to show that the 
energy inequality (|5]) holds in each it follows from (156|) and (l58|) that 


2p 


nqi 

2 


in 


Recall that p = pi and g = g^ in When B is antisymmetric, the energy inequality 
([5]) is reduced to 


= dt {x{pi) + + V ■ I {x{pi) + ^ +p(p 0) < 0 in x [0,T). 
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Direct computations yield 


= - ^dtp{pi) + dtlipi)^ + V ■ (lipi) + (1 - — 


n 


xm ■ 


m 


pi 


n 


■“Jo ,.pidtpi + m'’-\/pi b 

= IL \ dtxit) - x{t) - 2 -+ m ■ V 

^ I Pi 


2X{pi) + (2 - n)p{pi) 
npi 


2X{pi) + (2 - n)p{pi) , 

npi J J 

where V ■ = —dtPi in jg used in the last equality. Furthermore, the construction 

gives 

\dtpi\ = WT{t){po - P^i)\ < I^tWId||Vpo||l-(^(0) < CT-^P, 

|Vpi| < ||Vpo||L-(^(i)) < ^, P<Pi<P, |m^| = ^JnpiQi < 

therefore, the energy inequality fl85ll follows from 

xHt) < -CXQX^/^ - c^x - C{Qx"'^ - Co, 0 < t < T, (86) 


where the positive constants Cj = Cj{p, p) are independent of and the choices of the 
constants T and P. Therefore, choosing the constants T and P to be sufficiently small, 
we can hnd a smooth function y satisfying (1771) . (jUj), and fl86l) on [0, cxo). It is not hard 
to see that T can be chosen to be of the order of C{p,p)g~^. This hnishes the proof of 
the theorem. □ 
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